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Abstract 

Background: Genomic regions controlling abdominal fatness (AF) were studied in the Northeast Agricultural 
University broiler line divergently selected for AF. In this study, the chicken 60KSNP chip and extended haplotype 
homozygosity (EHH) test were used to detect genome-wide signatures of AF. 

Results: A total of 5357 and 5593 core regions were detected in the lean and fat lines, and 51 and 57 reached a 
significant level (P<0.01), respectively. A number of genes in the significant core regions, including RBI, BBS7, MAOA, 
MAOB, EHBPl, LRP2BP, LRPIB, MY07A, MY09A and PRPSAPl, were detected. These genes may be important for AF 
deposition in chickens. 

Conclusions: We provide a genome-wide map of selection signatures in the chicken genome, and make a 
contribution to the better understanding the mechanisms of selection for AF content in chickens. The selection for 
low AF in commercial breeding using this information will accelerate the breeding progress. 
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Background 

The linkage disequilibrium (LD) is important in livestock 
genetics for its key role in genomic selection [1] and 
detecting the causal mutations of economically import- 
ant traits [2-6]. Based on the LD information, there are 
two main methods to detect genes underlying pheno- 
typic variation, including one from phenotype to genome 
and another one from genome to phenotype. The first 
method is performed by targeting particular candidate 
genes or by quantitative trait loci (QTL) mapping and 
positional cloning of QTL. In the second method, pat- 
terns of LD in populations that are incompatible with 
the hypothesis of genetic neutrality are identified, and 
these patterns are selection signatures [7]. The aim of 
the second method is to identify artificial selections by 
statistically evaluating the genomic data [7] . 
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Allele frequencies underlying selection are expected to 
change. A neutral mutation will take many generations 
until the mutated allele reaches a high or low population 
frequency. In this case, the LD between the mutation 
and its neighboring loci will be degraded because of the 
recombination in every generation [8]. The frequency of 
a novel mutation will increase or decrease more rapidly 
than the neutral mutation because it is underlying artifi- 
cial selection, so that the surrounding conserved haplo- 
type was long [9,10]. This is the background of the 
extended haplotype homozygosity (EHH) statistic 
method used to detect selection signatures [11]. There 
are also many other methods to detect selective sweeps 
from DNA sequence data, including the Tajima's [12] 
and Fay and Wu's //-test [13] for selected mutations, 
measuring large allele-frequency differences among 
populations by Fst [14], and the integrated Haplotype 
Score (iHS) [15], which is an extension of the EHH stat- 
istic [11]. Among these methods, the EHH test is par- 
ticularly useful [7,11]. The EHH test is used to detect 
artificial selections according to the characteristics of 
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haplotypes within a single population, and do not re- 
quire the genotype of the ancestor [7]. Furthermore, the 
EHH test is less sensitive to ascertainment bias than 
other approaches, so it was designed to work with SNP 
rather than sequencing data [7,16]. 

The broilers used in this study were selected for 
eleven generations and genomic regions controlling AF 
deposition are expected to exhibit signatures of select- 
ive sweep. The aim of this study was to identify the se- 
lection signatures underlying the artificial selection for 
AF in chicken and to investigate the genes important 
for AF deposition. 

Methods 

Ethics statement 

All animal work was conducted according to the guide- 
lines for the care and use of experimental animals estab- 
lished by the Ministry of Science and Technology of the 
People's Republic of China (Approval number: 2006- 
398) and approved by the Laboratory Animal Manage- 
ment Committee of Northeast Agricultural University. 

DNA samples and data preparation 

Broilers used in this study were from two Northeast 
Agricultural University broiler lines divergently selected 
for AF content (NEAUHLF). The two lines have been 
selected since 1996 using AF percentage (%AFW or 
AFP) and plasma very low-density lipoprotein (VLDL) 
concentration as selection criteria [17]. The two lines 
were selected for 11 generations and the AFP changes 
over the 11 generations are shown in Figure 1. A total of 
475 individuals from generation 11 of NEAUHLF were 
used in this study. 

Genotyping was carried out using the Illumina 
chicken 60K SNP chip containing a total of 57636 
SNPs. Markers were filtered to exclude loci with un- 
known positions, monomorphic loci and loci with a 
minor allele frequency <0.05. 
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Figure 1 The separation of AFP over 1 1 generations between 
lean and fat lines. 



The haplotype and LD analysis 

The fastPHASE [18] (http://depts.washington.edu/fphase/ 
download/) was used to reconstruct the haplotypes for 
every chromosome using the default parameters. The 
reconstructed haplotypes were inserted into HAPLOVIEW 
v4.1 [19] to estimate LD statistics based on pairwise r^ and 
to construct the blocking pattern in the candidate regions 
of interest to enable selection signature analysis. 

The EHH test 

The "core region" was defined as the region in the gen- 
ome characterized by the strong LD among SNPs involv- 
ing a set of "core haplotypes" [7]. The Sweep v.1.1 (http:// 
www.soft82.com/get/ download/windows/sweep/) was 
used to identify the core regions [11]. The algorithm 
defined a pair of SNPs to be in strong LD if the upper 
95% confidence bound of D' is between 0.70 and 0.98 [20]. 
The program was set to select core regions with at least 
two SNPs. EHH was defined as the probability that two 
randomly chosen haplotypes carrying the candidate core 
haplotype were homozygous for the entire interval span- 
ning the core region to a given locus [11]. The EHH 
test [11] was based on one of the core haplotype vs. other 
haplotypes in the same position. The "Relative Extended 
Haplotype Homozygosity" (REHH) statistic corrects EHH 
for the variability in recombination rates [7]. It was com- 
puted by EHHt / EHH; with EHH defined as the decay of 
EHH on all other core haplotypes combined. The REHH 
value was used in the current study to determine the se- 
lection signatures. To determine the significance of REHH 
values, the haplotypes were ordered into 20 bins according 
to their frequencies [7]. The REHH values of each haplo- 
type in a candidate region were compared with all equally 
frequent haplotypes and the P-values were obtained [11]. 
The significant selection signatures were defined as 
P<0.01. 

Results 

Markers and core haplotypes 

A total of 43034 SNPs on 28 autosomes in chickens 
were included in the selection signature analysis 
(Table 1). These markers covered 950.68 Mb of the gen- 
ome, with an average of 22.09 kb between adjacent 
markers. 

For the SNPs analyzed in this study, the average minor 
allele frequency was 0.29 ± 0.13. A summary of genome- 
wide markers and core haplotype distribution in the data 
set is shown in Table 1. A total of 5357 and 5593 core 
regions spanning 549523.91 kb and 480784.79 kb of the 
genome, respectively, in the lean and fat lines were 
detected (Table 1). Mean core region length was esti- 
mated as 102.58±37.24 kb and 85.96±26.65 kb, with a 
maximum of 2288.64 kb and 2191.34 kb in the lean and 



Table 1 Summary of genome-wide marker and core region (CR) distribution in the lean and fat lines 
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^The number of SNPs. 

^Total length covered by core regions. 

^The proportion of total core region lengths on chromosome length. 
'^Number of SNPs In core regions. 

^The proportion of total number of SNPs in core regions on number of SNPs used. 
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Figure 2 Distribution of SNP numbers in the core regions (A) 
and the length of core regions (B) in lean and fat lines. 



fat lines, respectively (Table 1). Chromosome 1 was the 
largest chromosome in chickens, and it had the largest 
haplotypic structures in the genome, which covered 
110644.43 kb and 105728.03 kb in the lean and fat lines, 
respectively. For each chromosome, the proportion of 
length covered by core regions vs. total length, as well as 
the number of SNPs forming core regions vs. the total 
number of SNPs, are shown in Table 1. The distribution 
of the size of core regions is shown in Figure 2. Overall, 
25069 and 22180 SNPs in the lean and fat lines. 



respectively, participated in forming core regions, with a 
range of 2 to 19 SNPs per core. 



Whole genome selection signatures 

For all 5357 and 5593 core regions in the lean and fat 
lines, respectively, a total of 44822 and 46775 EHH tests, 
with an average of 8.37 and 8.36 tests per core region, 
were calculated. To find outlying core haplotypes, we 
calculated REHH at 1 Mb distances both on the up- 
stream and downstream sides. Figure 3 shows the distri- 
bution of REHH values vs. haplotype frequencies in the 
lean and fat lines, respectively. Corresponding P-values 
are indicated by different colored symbols. The -logio of 
the P-values associated with REHH against the chromo- 
somal position was plotted to visualize the chromosomal 
distribution of outlying core haplotypes with frequency 
<25% (Figure 4). The results indicated that these selec- 
tion signals were not uniformly distributed across all 
chromosomes, with a substantial overrepresentation on 
chromosomes 1, 2, 3 and 4. 

The genome-wide statistics of the selection signature 
test, including the number of tests and outlying core 
haplotypes for each chromosome, are presented in 
Table 2. Of 16677 and 18346 tests on core haplotypes 
with frequency >0.25, there were 51 and 57 tests with 
P<0.01 in the lean and fat lines, respectively. There were 
153 and 251 tests with P<0.05 in the lean and fat lines, 
respectively. 

The conformity of the distribution of Tukey's outliers 
was examined, with outlying core haplotypes defined at 
the threshold level of 0.01. Figure 5 displays box plots of 
the distribution of -loglO (P-values) within each bin of 
core haplotype frequency. The results indicated that the 
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Figure 3 Distribution of REHH vs. core haplotype frequencies in the lean and fat lines. Core haplotypes with P-values lower than 0 05 and 
0.01 are presented in blue and red, respectively. 
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Figure 4 Genome-wide map of P-values for core haplotypes with frequency >0.25 in lean and fat lines, respectively. Dashed lines 
display the threshold level of 0.01 . 
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Table 2 The number of tests on core haplotypes (CH) (both sides) with frequency>0.25 and P-values of REHH test 
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Figure 5 Box plot of the distribution of P-values in core haplotype frequency bins in the lean (left) and fat (right) lines. The dashed and 
continuous lines indicated the threshold P-values of 0.01 and 0.001, respectively. 
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Table 3 Statistics summary for core haplotypes with P<0.01 after the relative extended haplotype homozygosity 
(REHH) test 
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y^ZUUIIZO, L-n/VAL/ 


3 


1 09^709/^-1 n4^4Q(nQ 
1 UZ3 / zsj \ UH-DH-yoy 


0 55 


049 


2 49 


0 001 9 


r t / vn/ vo 




1 4RQ^9Qn-1 4Q^7rm7 
1 '+oyjzyu 1 ^+^3 / U3 / 


0 44 


0 72 


3 1 6 


n 004R 

U.UU'+O 


PLCB4 


3 


9ftQS7^4Q-9(^QQfnfn1 R 


045 


080 


3 69 


0 0009 


1 


3 


9ftQS7S4Q-9f^QQf^f^1 R 


045 


083 


3 47 


0 001 3 


j 


3 


97'^n'^Rnn-97'^'^^m n 

Z/ 3U30UU Z/ 3333 1 U 


0 52 


0 92 


3 78 


n oooQ 
u.uuuy 


SRBD1 


3 


9 7'^R9QQ^-9 74^00(^7 
Z / DOZ.^yIi Z. 1 '+3UUO / 


0 54 


0 84 


3 77 


n OOOQ 
u.uuuy 


SRBDl 




^^^^^71 R-^^ft1 n4fnfn 
33333/ lO 33DIU'-KJU 


0 47 


0 66 


3 02 


0 003 1 


E 1 C4G2 


3 


ftRQ^ft^9n-f^Qn7f^99^ 

U0^3U3ZU / KJjLjLJ) 


0 27 


0 98 


3 79 


0 0041 


RPF? GTF^Cf) n^F4R4 CDKIQ 


4 


^S99^SQ-^SS14Q4 


0 59 


0 54 


3 1 0 


0 0023 


MRMI ^ 

IVIUI VI ) 


4 


Q'^fiR7fi 1 -Qf^n4R7 1 

y300 / U 1 yOUH-O / 1 


0 55 


1 00 


2 95 


0 0040 




4 


1 llfJ^f\QRA 7R1 Q^^4 


041 


1 00 


3 56 


0 0037 




4 


Af^^ 4Q1 1 ft-4(^1 Qn97Q 


0 35 


0 99 


3 35 


0 0044 


FRFr, DfiA^M^ 


4 


c; 4. 7 4.4. Q 0 - ^ ^ 4 7 9 9 n Q 

33t^Zt^t^OU 33t^/ ZZUy 


0 55 


0 9 


2 60 


0 0005 


( nr > 1 1 ooj/ 


4 


R'^O'^I f^'^7-R^ 1 1 7Q74 

03U3 1 U3 / 03 1 1 / y / 


0 39 


0 87 


3 72 


0 0022 


1 
1 


c: 
J 


Q74nQ41 -QR9R144 
y / H-uyt 1 70ZO 1 


0 49 


1 00 


3 42 


0 001 4 


IFACnl CTRQ MR\/n 




9'^R9^1 1 ^-9'^R791 R7 
Z30Z3 1 1 3 Z30/ Z 1 O/ 


041 


0 78 


3 1 9 


0 0030 


093582 




49^09^1 7-49(^7Q4(nn 
H-Z3yZ3 1 / H-ZD/ y^-ou 


0 30 


1 00 


3 70 


n 0044 


1 


5 


fn^4fnfnn 1 -f^f^9f=i 1 4S 

\J3t^UUiJ 1 UUZU 1 T^3 


0 39 


1 00 


442 


0 001 1 


1 


6 


^^^^44SQ-^S^Qn^4(^ 


0 38 


0 99 


3 72 


0 0030 


PTPRE 




19.9.fS^f,f\A-19.QC\f{\AA 
zoooyoo'H- zoyuu3H-H- 


0 31 


1 00 


4 88 


0 001 5 


MYLK 




■^'^f^74nQR-^^71 ^1 99 
33Q / H-Uyo 33 / 1 3 1 ZZ 


0 57 


048 


2 09 


0 0027 




8 


6107407-617910^ 


0.35 


0.63 


3.72 


0.0020 


IFR^ KIAA1614 XPRl 


9 


1 6264832-1 6366749 


0.45 


0.97 


4.04 


0.0024 


PSMD1,ARMC9, B3GNT7 


10 


5831963-5856349 


0.59 


0.97 


1.83 


0.0034 


/ 


10 


19717086-19745274 


0.48 


1.00 


2.90 


0.001 1 


CHSYl 


n 


17094961-17160195 


0.30 


0.63 


3.35 


0.0001 


BCDOl, GAN 


13 


2628777-2664596 


0.35 


1.00 


3.94 


0.0048 


Q5ZHQ9 


13 


2726706-2746894 


0.39 


0.79 


4.06 


0.0041 


1 


13 


16758621-16783127 


0.25 


0.69 


3.26 


0.0050 


F5TL4 


15 


7345639-7377799 


0.25 


0.54 


2.38 


0.0004 


SEZ6L, ASPHD2, HPS4 


18 


9949736-10015444 


0.39 


1.00 


2.34 


0.0031 


SPAG9 
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Table 3 Statistics summary for core haplotypes with P<0.01 after the relative extended haplotype homozygosity 
(REHH) test (Continued) 



18 


10117401-10135964 


0.39 


1.00 


2.64 


0.0013 


F1NM51 


19 


8727596-8786448 


0.60 


0.55 


1.77 


0.0038 


Msn 


20 


9090808-9113453 


0.29 


0.95 


6.18 


0.0036 


MYTl 


20 


9246245-9278998 


0.32 


0.80 


5.39 


0.0040 


E1C8M0 


20 


9879361-9899719 


0.27 


0.96 


6.40 


0.0030 


CSK21 


22 


2952274-3002268 


0.29 


0.29 


3.93 


0.0039 


/ 


Fat line 






0 33 


0 81 


445 


0 001 8 


TCF20 




JO 1 zuuuy joz 1 D joh- 


0 55 


0 99 


3 1 5 


0 001 6 


DRi i\/nd 




OU 1 / 1 U / O DUZ J'4-/ / 1 


0 46 


1 00 


4 1 2 




1 

1 




O/ / OjOOZ D/ O jUO I O 


0 26 


0 96 


4 06 


0 0026 


1 




QOZ 1 JQ 1 / OOZ J / ZH 1 


0 61 


0 94 


4 21 


0 001 6 


S0X5 




f\Qf^'?A^ Rft-(^Q(^S(^'^c;7 

0^03^+ 1 OO DyOOU J J / 


0 66 


0 99 


2 89 


0 0003 


1 





1 ni S^Sf=i1 S-inifn^Sfnfn7 


0 29 


1 00 


5 24 


0 0005 


C/1A/JCM/ 

3A1/ VI3/ V / 




1 147RQAR7-1 14.R7^(=^9^ 


0 29 


0 99 


3 66 


n 004R 


IViMKJD, IVir\\Jf\ 




1 1cjQOQQQ^-l 7ftni 1 QR4. 


0 35 


1 00 


3 77 


0 0036 


E 1 BTB5 




\ Jt-OOJ J 1 U 1 D'-t/ JAy\jJ 


0 72 


0 89 


2 09 


0 0034 






1 R1 Rnn997-1 R1 RR^^4^ 

1 O 1 OUUZZ / 1 O 1 OO J DH-J 


0 33 


0 99 


4 00 


0 0033 


r\ 1 A\JVU 





1 R1 Rnn997-1 R1 RR^S4S 

1 O 1 OUUZZ / 1 O 1 OO J Jt^ J 


0 33 


1 00 


3 81 


0 0043 




2 


7ft7ftRR41 -7ftR^4^9'^ 
/ U / DoOH- 1 / Do DH-JZ J 


0 3 1 


1 00 


4 30 


0 0041 




2 


i^i9n^Q'^^-i^i9mrmQ 

1 J 1 zujyj J 1 J 1 ZD 1 ujy 


0 62 


0 82 


2 32 


0 0033 


TRAPPC9 


2 


1 ^■^l 1 7009-1 ^■^l zl'^RR'^ 
1 J 3 1 1 / uyz 1 J J 1 H-joo J 


0 71 


0 77 


1 68 


0 0038 




3 


Q1 77Qn7-Q799R9^ 


0 38 


0 95 


4 00 


0 0024 


EHBP1 




Q1 77Qn7-Q799R9^ 


0 38 


0 90 


3 66 


0 0039 


EHBP1 


3 
0 




0 30 


0 98 


4 98 


0 001 6 


1 


3 


94Q4^R'^Q-94QR(^779 
ZH-yH-jojy ZH-yoo/ / z 


0 61 


0 70 


2 1 5 


0 001 4 


1 
1 


3 


AAlf^^ 1 1 f^-AA'^ 1 1 4.Q'^ 
H-H-zoD 1 1 o H-H-j 1 1 H-y J 


0 40 


0 99 


3 50 


0 0050 




3 


fiQRfi'^R^n-f^QQnfnfnQR 
oyoojoDU oyyuooyo 


0 34 


1 00 


4 29 


0 0038 


1 


3 


R^R741 ^7-R^Q^14.7^ 


041 


1 00 


3 97 


0 0031 




3 


Q7997f^Rn-Q7^^7Qnfn 
y / zz / Qou y / 33/ yuo 


0 28 


0 99 


3 81 


0 0042 


1 
1 


4 


1 1 ^R9 1 4.1 -1 1 i^49c;3R 

1 1 JOZ 1 H 1 1 1 Q'+ZJjO 


0 27 


0 96 


4 53 


0 0029 


1 


4 


4nft^^^Q^-4n71 ^404 

T^UU333yj T^U / 1 3t^Ut^ 


0 32 


0 93 


440 


0 001 0 


Cdnrf 7(1 1 RP7RP ^MY?^ 


4 


55950677-55991394 


0.28 


0.99 


4.55 


0.0028 


1 


4 


55950677-55991394 


0.28 


1.00 


4.36 


0.0036 


1 


4 


86719441-86754976 


0.57 


0.75 


2.88 


0.0048 


1 


5 


556571-628531 


0.25 


0.67 


4.37 


0.0026 


nNYX6, PLCB2, BUBW, PAK6 


5 


40239840-40261525 


0.29 


0.97 


4.48 


0.0038 


V5X2, F1N9P5 


5 


47240577-47282933 


0.40 


0.95 


2.91 


0.0029 


RIN3, LQMN 


5 


59811459-59880511 


0.39 


0.66 


3.36 


0.0041 


1 


6 


26756202-26793956 


0.34 


0.98 


3.94 


0.0049 


1 


6 


29341938-29401207 


0.32 


0.97 


5.50 


0.0007 


ABLIMl 


7 


30090927-30155133 


0.30 


0.88 


3.77 


0.0033 


F1NF72 


7 


31374271-31418061 


0.42 


1.00 


2.88 


0.0049 


LYPDl, NCKAP5 


7 


33795201-33904515 


0.47 


1.00 


2.20 


0.0025 


LRP1B 


7 


36818722-36875768 


0.26 


0.99 


3.95 


0.0024 


Q9DEH4 


7 


37031922-37124566 


0.56 


1.00 


2.92 


0.0047 


STAM2, FMNL2 


8 


5597-492518 


0.56 


0.99 


1.94 


0.0036 


F INF 53 
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Table 3 Statistics summary for core haplotypes with P<0.01 after the relative extended haplotype homozygosity 
(REHH) test (Continued) 



g 




047 


0 92 


4 20 


0 0004 


NEK7 


9 


!7 J Z. !7 1 J UU / V J 


041 


0 71 


3 85 


0 0044 


/ 


1 0 




0 50 


0 78 


5 04 


n 0004 
u.uuu^ 


MYDQA FWDMA 


1 1 




0 47 


0 52 


3 1 1 


0 001 4 


/ 

I 


1 1 




0 58 


0 35 


1 80 


0 0040 


1 


1 2 


1 1 c;71 QQ-I 1 7nifnQ 


0 37 


063 


4 57 


0 0030 


1 


1 3 


1 J J J J Ji. 1 Ut^W 1 Jt^ 


040 


1 00 


3 41 


0 0049 




14 


Rn4fin^Q-R1 7^fn9Q 
out-ouj" o 1 / jozy 


0 42 


0 86 


3 09 


0 001 5 


1 

I 


15 


8495796-8543001 


029 


1.00 


4.86 


0.0003 


TBX6, CRKL, KLHL22 


15 


8495796-8543001 


0.29 


0.99 


4.43 


0.0006 


TBX6, CRKL, KLHL22 


17 


3250605-3271593 


0.27 


0.97 


3.91 


0.0033 


1 


17 


4062173-4087131 


0.26 


0.94 


3.77 


0.0040 


C4PCF3 


18 


4433126-4445816 


0.26 


0.79 


5.99 


0.0037 


PRPSAP1 


18 


8365846-8400245 


0.47 


0.80 


3.53 


0.0038 


1 


23 


935267-970086 


0.31 


0.91 


5.40 


0.0040 


EDN2 


24 


5613517-5633477 


0.28 


0.86 


4.38 


0.0024 


ZWIO, FINCIO 


24 


6145308-6158962 


0.31 


0.80 


5.30 


0.0047 


/ 



^REHH and P-values are presented for upstream and downstream sides from each core haplotype, respectively. 



extreme outliers appear in the small haplotype frequen- 
cies bins. 

Mapping selection signatures to genes 

A summary of statistics for 51 and 57 positively selected 
core regions with P<0.01 of the REHH tests in the lean 
and fat lines, respectively, is presented in Table 3. Corre- 
sponding genes were identified by aligning the core posi- 
tions with the chicken genome sequence (Table 3). The 
full genes names were from Ensembl online (http:// 
www.ensembl.org/index.html). A total of 66 and 46 
genes in the core regions were detected in the lean and 
fat lines, respectively, including RBI (retinoblastoma 1), 
BBS7 (Bardet-Biedl syndrome 7), MAOA (monoamine 
oxidase A), MAOB (monoamine oxidase B), EHBPl (EH 
domain binding protein 1), LRP2BP (LRP2 binding pro- 
tein), LRPIB (low-density lipoprotein receptor-related 
protein IB), MY07A (myosin VIIA), MY09A (myosin 
IXA) and PRPSAPl (phosphoribosyl pyrophosphate 
synthetase-associated protein 1). The haplotype analysis 
of these genes revealed that the haplotype frequencies 
were significantly different (P<0.01) between the two 
lines (Table 4). 

Mapping selection signatures to QTLs 

The chicken QTL database available online (http:// 
www.animalgenome.org/ cgi-bin/QTLdb/GG/index) was 
explored to identify any overlapping of the core 
regions with significant REHH P-values (P<0.01) and 
published QTLs in chickens. The approximate posi- 
tions of the overlapping QTLs for each core region are 



listed in Table 5. There were many overlaps between 
the core regions with significant REHH P-values 
(P<0.01) and published QTLs for AF content in 
chickens. 

Discussion 

Selective sweep is used to detect genomic regions with 
reduced variation in allele frequency in any population 
experiencing divergent selection for specific traits. Here, 
we determined the feasibility of the selective sweep ap- 
proach for finding genes important for AF deposition in 
chickens. The long-range haplotype test was employed, 
which detects selection signature by measuring the char- 
acteristics of haplotypes within the lean and fat lines di- 
vergently selected for AF content. There were 5357 and 
5593 core regions in the lean and fat lines, respectively. 
When comparing the average marker spacing with mean 
core length and number of SNPs forming cores, we 
revealed that core regions are more likely to appear in 
regions with higher marker density. 

The selection signatures on the whole genome were 
calculated, and a subset of putative core regions with 
significant REHH P-values (P<0.01) was identified. The 
genes in these core regions were detected and 10 genes, 
including RBI, BBS7, MAOA, MAOB, EHBPl, LRP2BP 
LRPIB, MY07A, MY09A and PRPSAPl, were important 
for fatness. Among these 10 important genes, seven 
genes, including RBI, BBS7, MAOA, MAOB, EHBPl, 
LRP2BP and LRPIB, were all in the QTL regions 
reported previously for AF in chickens (Table 5). Al- 
though the other three genes, including MY07A, 
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Table 4 Haplotype frequencies in the lean and fat lines of the core regions including 10 important genes 



Gene and core regions 


Haplotype Number 


Haploptypes 


Haplotype frequency 
Lean line Fat 


line 


n 1 1 

P-value 


MAOB, MAO A 


1 


CAAGG 


0.645 


0.615 


<0.001 


Chrl: 1 14789487-1 14875623 




AAAGA 


0.197 


0 






3 


CGGAG 


0.158 


0.269 






4 


CGAGA 


0 


0.077 






5 


AAAGG 


0 


0.038 




RBI 


1 


GGAA 


0.421 


0.410 


<0.001 


Chrl: 173098805-173190831 


^ 


GAGG 


0.368 


0.103 






3 


GAAA 


0.21 1 


0.192 






4 


AAGG 


0 


0.244 






5 


AGGA 


0 


0.038 






6 


GAGA 


0 


0.013 




MY07A 


1 


AGG 


0.618 


0.090 


<0.001 


Chrl: 198071099-198113519 


2 


GAA 


0.316 


0.207 






3 


GGA 


0.066 


0.652 






4 


GAG 


0 


0.037 






5 


GGG 


0 


0.014 




EHBP1 


1 


GGG 


0.855 


0.090 


<0.001 


Chr3: 91 77907-9222825 


2 


GAG 


0.132 


0.359 






3 


AGG 


0.01 3 


0.128 






■1 


GGA 


0 


0.423 




LRP2BP 


1 


GGGG 


0.443 


0.487 


<0.001 


Chr4: 40653593-40713404 


2 


AAAA 


0.338 


0.21 1 






3 


GGAA 


0.1 76 


0 






4 


AAGG 


0.044 


0.303 




BBS? 


1 


AGGC 


0.605 


0.282 


<0.001 


Chr4: 55424480-55472209 


- 


GAAA 


0.368 


0.301 






3 


AAAA 


0.026 


0 






4 


AGAC 


0 


0.198 






5 


AGAA 


0 


0.161 






6 


GAAC 


0 


0.058 




LRPIB 


1 


AGAGAC 


0.361 


0.013 


<0.001 


Chr7: 33795201-33904515 


2 


GGAGGA 


0.197 








3 


AGAAGA 


0.105 


0.154 






4 


GGGGGA 


0.066 


0.449 






5 


AGAAGC 


0.057 


0.346 






6 


GAGGGA 


0.055 


0.038 






7 


GAGAGA 


0.050 


0 






8 


GGAAGA 


0.049 


0 






9 


GAGGAA 


0.026 


0 






10 


GGAGAA 


0.018 


0 






11 


AGAAAC 


0.016 


0 




MY09A 


1 


GGGAA 


0.355 


0.051 


<0.001 


ChrlO: 763998-831991 


2 


AAGAA 


0.276 


0.358 






3 


AAGAG 


0.237 


0 






4 


GGGGA 


0.118 


0.013 






5 


AGGAA 


0.013 


0.065 
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Table 4 Haplotype frequencies in the lean and fat lines of the core regions including 10 important genes (Continued) 



6 AGAAA 0 0.500 





7 


AAAAA 


0 


0.013 




PRPSAP1 


1 


AGA 


0.816 


0.615 


<0.001 


Chr18: 4433126-4445816 


2 


GGG 


0.118 


0.026 






3 


AAG 


0.066 


0.090 






'1 


AGG 


0 


0.269 




^P-values of Fisher's Exact Test for difference 


analysis of haplotype frequencies 


between lean and fat lines. 









MY09A and PRPSAPl, were not in the QTL regions, 
these genes were also important for the AF deposition. 

The known functions of these 10 genes were analyzed 
and the results indicated that they were likely to be 
linked with fatness. The RBI gene regulates the C/EBP- 
DNA-binding activity during 3T3-L1 adipogenesis and 
plays a key role in adipocyte differentiation [40,41]. 

The BBS7 gene is a member of the Bardet-Biedl syn- 
drome (BBS) family. BBS is a pleiotropic genetic disorder 
characterized by obesity, photoreceptor degeneration, 
Polydactyly, hypogenitalism, renal abnormalities, and de- 
velopmental delay [42]. BBS is recognized to be a genet- 
ically heterogeneous autosomal recessive disorder mapped 
to eight loci [42]. Positional cloning and candidate genes 
identified six BBS genes, including BBSl, BBS2, BBS4, 
BBS6, BBS7, and BBSS [42]. These BBS genes may be im- 
portant for obesity. 

The MAO A and MAOB are two enzymes important 
for dopamine production. The dopamine levels influence 
the risk of obesity and MAOA and MOAB may be impli- 
cated in human obesity [43]. 

The EHBPl gene is required for insulin-stimulated 
GLUT4 movements [44]. Insulin stimulates glucose trans- 
port in adipose tissues by recruiting intracellular mem- 
brane vesicles containing the glucose transporter GLUT4 
to the plasma membrane [44]. The mechanisms involved 
in the biogenesis of these vesicles and their translocation 
to the cell surface were studied and the results indicated 
that EHDl and EHBPl are required for perinuclear 
localization of GLUT4, and the loss of EHBPl dis- 
rupts insulin-regulated GLUT4 recycling in cultured 
adipocytes [44]. This indicates that the EHBPl gene 
may be important in adipocyte differentiation. 

The LRP2BP and LRPIB genes are two members of 
the low-density lipoprotein receptor family that partici- 
pates in a wide range of physiological processes, includ- 
ing the regulation of lipid metabolism, protection 
against atherosclerosis, neurodevelopment, and transport 
of nutrients and vitamins [45] . 

The MY07A and MY09A are two myosin genes. A 
spontaneous mutant mouse line, Myo7'^'^^^'^\ was used 
to study the function of the MY07A gene, and the result 
indicated that the mutant male homozygous mice 
displayed decreased body weight and body fat [46]. 



The MY09A gene was in the BBS4 region of chromo- 
some 15q22-q23 [47], which might be important for 
obesity. 

The PRPSAPl gene is named as phosphoribosyl pyro- 
phosphate synthetase-associated protein 1. The results 
of differentially expressed genes associated with insulin 
resistance indicate that PRPSAPl gene is associated with 
percentage of body fat [48]. 

The associations of these 10 genes with obesity or lipid 
metabolism were mainly in humans and mice. Because 
of the high conservation of these genes between 
humans, mice and chickens, the 10 genes might also be 
important for AF deposition in chickens. 

Conclusions 

Our results provide a genome-wide map of selection sig- 
natures in two chicken lines divergently selected for AF 
content. There were 51 and 57 core regions showing sig- 
nificant P-values (P<0.01) of selection signatures in the 
lean and fat lines, respectively. In these core regions 
there were a number of important genes, including RBI, 
BBS7, MAOA, MAOB, EHBPl, LRP2BP LRPIB, 
MY07A, MY09A and PRPSAPl. These genes are im- 
portant for AF deposition in chickens. 
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Table 5 Reported QTL near the core regions with P <0.01 in the lean and fat lines 


Lean line 


Chr Core region (bp) 


Trait 


QTL position (bp) 


F-ratio 


P-value 


Reference 


1 39360501-39455853 


AFP 


1 937738-52700434 


1.474 


Suggestive 


[21] 


1 49926970-49964278 


AFP 


25998723-65961966 


1.732 


Suggestive 


[21] 




AFW 


25998723-65961966 


1.882 


Suggestive 


[21] 




AFW 


48175152- 51977642 


8.14 


Significant 


[22] 


1 173098805-173190831 


AFW 


158352237- 182910620 


3.18 


Significant 


[23] 




AFP 


171224834- 174526878 


20.34 


Significant 


[23] 


2 3631683-3739002 


AFW 


3097660- 4097660 


3.38 


Suggestive 


[24] 


3 3794973-3861882 


AFP 


800029-11 0574691 


1.364 


Suggestive 


[21] 


3 10257926-10454969 


AFP 


800029-1 10574691 


1.364 


Suggestive 


[21] 




AFW 


6841859-13986734 


8.16 


Significant 


[22] 




AFP 


6841859- 57396057 


7.9 


Significant 


[25] 




AFW 


6841859- 44850897 


74 


Significant 


[25] 


3 14895290-14957057 


AFP 


800029-1 10574691 


1.364 


Suggestive 


[21] 




AFW 


6841859-13986734 


8.16 


Significant 


[22] 




AFP 


6841859- 57396057 


7.9 


Significant 


[25] 




AFW 


6841859- 44850897 


74 


Significant 


[25] 




AFW 


1 3986734-25508863 


\ 


Suggestive 


[26] 


3 26957549-26996618 


AFP 


800029-1 10574691 


1.364 


Suggestive 


[21] 




AFP 


6841859- 57396057 


7.9 


Significant 


[25] 




AFW 


6841859- 44850897 


7.4 


Significant 


[25] 




AFW 


24160710-51592221 


\ 


Suggestive 


[27] 




AFW 


25508863- 35512024 


\ 


Suggestive 


[26] 


3 27303800-27335510 


AFP 


800029-1 10574691 


1.364 


Suggestive 


[21] 




AFP 


6841859- 57396057 


7.9 


Significant 


[25] 




AFW 


6841859- 44850897 


7.4 


Significant 


[25] 




AFW 


24160710-51592221 


\ 


Suggestive 


[27] 




AFW 


25508863- 35512024 


\ 


Suggestive 


[26] 


3 27382993-27430067 


AFP 


800029-1 10574691 


1.364 


Suggestive 


[21] 




AFP 


6841859- 57396057 


7.9 


Significant 


[25] 




AFW 


6841859- 44850897 


74 


Significant 


[25] 




AFW 


24160710-51592221 


\ 


Suggestive 


[27] 




AFW 


25508863- 35512024 


\ 


Suggestive 


[26] 


3 35555718-35610466 


AFP 


800029-1 10574691 


1.364 


Suggestive 


[21] 




AFW 


35512024-40755790 


18.5 


Significant 


[28] 




AFP 


35512024-40755790 


13.1 


Significant 


[28] 


4 17765695-17819334 


AFW 


17425871-18425871 


\ 


Significant 


[29] 


4 46149116-46190279 


AFW 


42005559- 51609571 


2.26 


Suggestive 


[30] 


4 55424480-55472209 


AFP 


51266614- 88408499 


16.0 


Significant 


[25] 


4 83051637-83117974 


AFP 


51266614- 88408499 


16.0 


Significant 


[25] 




AFW 


80258156-88408499 


6.9 


Significant 


[25] 




AFW 


81539616-84618310 


2.04 


Suggestive 


[30] 


5 23825115-23872187 


AFW 


18412554-42717839 


21.8 


Significant 


[25] 




AFP 


18723157- 43339045 


194 


Significant 


[25] 




AFW 


19782191- 30162990 


\ 


Suggestive 


[26] 




AFW 


19782191- 30162990 


7.04 


Significant 


[31] 
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Table 5 Reported QTL near the core regions with P <0.01 in the lean and fat lines (Continued) 


5 42592517-42679460 


AFW 


18412554-42717839 


21.8 


Significant 


[25] 




AFP 


18723157- 43339045 


19.4 


Significant 


[25] 




AFW 


37226264-53779276 


6.74 


Significant 


[31] 


6 35354459-35390346 


AFP 


29647151- 37399694 


6.9 


Significant 


[25] 


7 28869664-28906344 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 




AFW 


28166221- 29166221 


9.78 


Significant 


[32] 




AFW 


28166221- 29166221 


\ 


Significant 


[33] 


7 35674098-35715122 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 


9 16264832-16366749 


AFW 


13658592-23770679 


5.03 


Suggestive 


[22] 




AFW 


15457880-16457880 


7.0 


Suggestive 


[34] 


10 19717086-19745274 


AFP 


16519830- 20778533 


9.9 


Significant 


[28] 


13 16758621-16783127 


AFW 


16327806- 18173123 


2.10 


Suggestive 


[30] 


15 7345639-7377799 


AFW 


1917251- 10769106 


10.2 


Significant 


[25] 




AFP 


2388961-10769106 


12.8 


Significant 


[25] 




AFW 


2798507-10769106 


8.13 


Significant 


[22] 




AFW 


2798507-10769106 


5.67 


Suggestive 


[22] 




AFW 


3717446-7928397 


2.21 


Suggestive 


[30] 




AFP 


3717446-7928397 


2.22 


Suggestive 


[30] 


Fat line 


1 51248496-51279543 


AFP 


1 937738-52700434 


1.474 


Suggestive 


[21] 




AFW 


48175152- 51977642 


8.14 


Significant 


[22] 




AFP 


25998723- 65961966 


1.732 


Suggestive 


[21] 




AFW 


25998723- 65961966 


1.882 


Suggestive 


[21] 


1 58120009-58215364 


AFP 


25998723- 65961966 


1.732 


Suggestive 


[21] 




AFW 


25998723- 65961966 


1.882 


Suggestive 


[21] 




AFW 


55261695-67128747 


12.18 


Significant 


[35] 


1 60171076-60254771 


AFP 


25998723- 65961966 


1.732 


Suggestive 


[21] 




AFW 


25998723- 65961966 


1.882 


Suggestive 


[21] 




AFW 


55261695-67128747 


12.18 


Significant 


[35] 


1 67763862-67830818 


AFW 


67327367-68327367 


\ 


Significant 


[33] 


1 68213617-68257241 


AFW 


67327367-68327367 


\ 


Significant 


[33] 


1 101535615-101635667 


AFW 


89938943-167462479 


94 


Significant 


[36] 




AFW 


94157976- 102460326 


6.11 


Suggestive 


[35] 


1 114789487-114875623 


AFW 


113344161- 132660888 


7.90 


Suggestive 


[35] 




AFW 


114143603- 115143603 


7.1 


Significant 


[36] 


1 125909995-126011984 


AFW 


113344161- 132660888 


7.90 


Suggestive 


[35] 


1 181800227-181883545 


AFW 


158352237-182910620 


3.18 


Significant 


[23] 


3 9177907-9222825 


AFP 


800029- 110574691 


1.364 


Suggestive 


[21] 




AFW 


6841859- 13986734 


8.16 


Significant 


[22] 




AFW 


6841859- 13986734 


5.8 


Suggestive 


[22] 




AFP 


6841859-57396057 


7.9 


Significant 


[25] 




AFW 


6841859-44850897 


74 


Significant 


[25] 


3 16143474-16194865 


AFP 


800029- 110574691 


1.364 


Suggestive 


[21] 




AFP 


6841859-57396057 


7.9 


Significant 


[25] 




AFW 


6841859-44850897 


74 


Significant 


[25] 




AFW 


1 3986734-25508863 


\ 


Suggestive 


[26] 
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Table 5 Reported QTL near the core regions with P <0.01 in the lean and fat lines (Continued) 


3 24945839-24986772 


AFP 


800029- 1 10574691 


1.364 


Suggestive 


[21] 




AFP 


6841859-57396057 


7.9 


Significant 


[25] 




AFW 


6841859-44850897 


7.4 


Significant 


[25] 




AFW 


1 3986734-25508863 


\ 


Suggestive 


[26] 




AFW 


24160710-51592221 


\ 


Suggestive 


[27] 


3 44265116-44311493 


AFP 


800029- 1 10574691 


1.364 


Suggestive 


[21] 




AFP 


6841859-57396057 


7.9 


Significant 


[25] 




AFW 


6841859-44850897 


7.4 


Significant 


[25] 




AFW 


24160710-51592221 


\ 


Suggestive 


[27] 




AFW 


40755790-45203763 


7.5 


Significant 


[28] 




AFP 


40755790-45203763 


10.8 


Significant 


[28] 


3 69863850-59906698 


AFP 


800029- 1 10574691 


1.364 


Suggestive 


[21] 


3 85874137-85931473 


AFP 


800029- 1 10574691 


1.364 


Suggestive 


[21] 


3 97227680-97337906 


AFP 


800029- 1 10574691 


1.364 


Suggestive 


[21] 


4 40653593-40713404 


AFP 


40473174-41473174 


\ 


Significant 


[32] 


4 55950677-55991394 


AFP 


51266614- 88408499 


16.0 


Significant 


[25] 


4 86719441-86754976 


AFP 


51266614- 88408499 


16.0 


Significant 


[25] 




AFW 


80258156-88408499 


6.9 


Significant 


[25] 


5 40239840-40261525 


AFW 


18412554- 42717839 


21.8 


Significant 


[25] 




AFP 


18723157- 43339045 


19.4 


Significant 


[25] 




AFW 


37226264- 53779276 


6.74 


Significant 


[31] 




AFW 


40158255- 41 158255 


\ 


Significant 


[37] 




AFW 


40158255- 41158255 




Significant 


[38] 




AFP 


40158255- 41 158255 


\ 


Significant 


[38] 


5 47240577-47282933 


AFW 


37226264- 53779276 


6.74 


Significant 


[31] 


5 59811459-59880511 


AFW 


51748760-60234891 


\ 


Significant 


[26] 




AFW 


53867807-62098509 


1 1.87 


Significant 


[31] 




AFW 


53867807-62098509 


6.82 


Significant 


[31] 


7 30090927-30155133 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 


7 31374271-31418061 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 


7 33795201-33904515 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 




AFW 


32440861-34526547 


2.08 


Suggestive 


[30] 


7 36818722-36875768 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 


7 37031922-37124566 


AFW 


25306930- 38010856 


\ 


Suggestive 


[27] 


9 2952291-3007034 


AFW 


2798942-3798942 


\ 


Significant 


[32] 




AFP 


2972071-3972071 


\ 


Significant 


[32] 


1 1 9804894-9826761 


A.FW 


6272742- 12810705 


2.15 


Suggestive 


[30] 


IT 11 c-71 nn 1 1 "7ni ^^n 

1 z 1 1 5 / 1 yy- 1 1 /u 1 6y 


AFP 


/34z0y- I22/5U26 


5.22 


Significant 


L28J 




AFP 


734209- 1 2275026 


4.51 


Significant 


[28] 




AFP 


813709-1813709 


\ 


Significant 


[32] 


15 8495796-8543001 


AFW 


1917251- 10769106 


10.2 


Significant 


[25] 




AFP 


2388961- 10769106 


12.8 


Significant 


[25] 




AFW 


2798507- 10769106 


8.13 


Significant 


[22] 




AFW 


2798507- 10769106 


5.67 


Suggestive 


[22] 


23 935267-970086 


AFW 


74802-1074802 


\ 


Significant 


[39] 
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